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ABSTRACT 

The variation of the kinematical properties of the Galactic thick disk with Galactic 
height Z are studied by means of 412 red giants observed in the direction of the south 
Galactic pole up to 4.5 kpc from the plane. We confirm the non-null mean radial motion 
toward the Galactic anticenter found by other authors, but we find that it changes sign 
at 1-2^1=3 kpc, and the proposed inward motion of the LSR alone cannot explain these 
observations. The rotational velocity decreases with \Z\ by —30 km s~^ kpc~^, but the 
data are better represented by a power-law with index 1.25, similar to that proposed 
from the analysis of SDSS data. All the velocity dispersions increase with |Z|, but the 
vertical gradients are small. The dispersions grow proportionally, with no significant 
variation of the anisotropy. The ratio (Tu/o"w=2 suggests that the thick disk could 
have formed from a low-latitude merging event. The vertex deviation increases with 
Galactic height, reaching ~20 degrees at |Z|=3.5 kpc. The tilt angle also increases, and 
the orientation of the ellipsoid in the radial-vertical plane is constantly intermediate 
between the alignment with the cylindrical and the spherical coordinate systems. The 
tilt angle at \Z\=2 kpc coincides with the expectations of MOND, but an extension of the 
calculations to higher \Z\ is required to perform a conclusive test. Finally, between 2.5 
and 3.5 kpc we detect deviations from the linear trend of many kinematical quantities, 
suggesting that some kinematical substructure could be present. 

Subject headings: Galaxy: disk - Galaxy: kinematics and dynamics - Galaxy: structure 
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1. INTRODUCTION 

The formation mechanism of the thick disk is one of the most noticeable grey points of om' 
current understanding of the Galactic formation and evolution process. This uncertainty is par- 
ticularly unfortunate because, being the thic k disk a ubiquitous feature among spiral galaxies 
(jPalcanton &: Bernsteinll2002l : ISeth et alJl2005l ). its formation must be a common s tage in the early 
evolu tion of late- type galaxies. During the nearly three decades since its discovery (iGilmore &: Reid 
19831 ). many models have been proposed. The merging scenario, in which the early merging of a 
small satellite galax y heats a primord i al Galactic disk pr oducing an old, thick, and kinematically 
hot disk population (iQuinn et al.lll993l : IWalker et al.lll996l ). has been the most favored model in the 



last decade, following the evidence that the thin and thick disks are chemically distinct (jFuhrmann 



1998; Gratton et al 



ample, 



2001 



Bournaud et al 



Bournaud et al 



Nevertheless, even this scenario is not free of problems (see, for ex- 



200911. and alternative inodels have recently been drawn to attent ion (e.g., 
2OO9I : ISchonrich fc Binn"evll2009l : iRoskar et aPbood : kssmann et~al]l201lh . 



It is clear that the general properties of the Galactic thick disk, such as its mean metallicity 
or mean kinematics, are not enough to distinguish between the models proposed for its formation. 
Moreover, the merging scenario has many free parameters, such as the mass of the merging satellite 
and the inclination of its orbit with respect to the Galactic disk, and the observations must constrain 
them if the quality of the model is to be finally assessed. In the last few years, theoretical simulations 
have started to cast predictions of the stellar distribution of stars in the multi-dime nsional spatial- 



kinematical-chemica l space (e.g.. iLoebman et al.ll201ll : lKobavashi &: Nakasatdl201ll ). For example 



Gomez et al.l (|201ll ) have shown that, within the merging scenario, the time of impact, and the 



mass and orbit of the infalling satellite c an be deduced frorn the distribu ti on of the expected kine - 



mat ical clumps in the U -V plane, while IVillalobos fc Helmil (120081 . l2009l ). IVillalobos et al.l ()2010l ) 



and iPurcell et al.l ([20091) find that the a\j / uw ratio and its variation w ith Galactocentric d i stance 



are excellent indicators of the inclination angle of the merging event. IVillalobos Helmil ()2008l ) 



and lBekki &: Tsujimotol (|201ll ) have also shown that the observed decrease of rotation velocity with 



distance from the plane points to a low-latitude merging. At the same time, the observations are 



tribution (e.g. Ilvezic et al 



2008 



Fuchs et al. 


2009: 


Bond et al. 


2010: 


Casetti-Dinescu et al. 


2011 



Detailed results of this kind are strongly needed because, through comparison with the expectations 
from the different models, they can be key to solve the puzzle of the Galactic thick disk formation. 

We are undertaking an extensive survey aimed to rev eal the kinematical and chemical vertical 
structure of the Galactic thick disk (ICarraro et al.ll2005l). P reliminarv kinematical results were 
presented by iMoni Bidin et al.l (j2008l ) and iMoni BidinI (|2009l ) , while the sample was also used to 



^Based on observations collected at the European Organization for Astronomical Research in the Southern Hemi- 
sphere, Chile (proposal IDs 075.B-0459(A), 077.B-0348(A)) 

^This paper includes data gathered with the 6.5-meter Magellan and the duPont Telescopes, located at Las 
Campanas Observatory, Chile. 
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investigate the Galactic dark d isk (jMoni Bidin et alj 120101 ) . and Lithium-rich stars in the Galactic 
thick disk (jMonaco et al.ll201ll ). In this paper, we focus on the kinematical results, studying the 
trend of kinematical quantities as a function of distance from the Galactic plane. In later papers of 
this series, the collected spectra will be used to measure the metallicity and chemical abundances 
of the sample, to study the variation of the thick disk chemistry with Galactic height. 



THE SAMPLE 



Our investigation is based on the sample of '^1200 red giants defined bv lGirard et al.l (|2006l ). 
vertically distributed with respect to the Galactic plane in a cone of fifteen degrees radius centered 
on the South Galactic Pole. The sample covers a large range of Galactic heights, from the plane to 
\ Z\ ~5 kpc, and it is volume-complete up to |.Z'|=3 kpc. All the objects have 2MASS photometr y 
(jSkrutskie et al.ll2006l ). and absolute proper motions from the SPM3 catalog (jGirard et al.l 12004 ). 
The sample was defined by the color cut 0.7< {J— Kg) <1.1, to isolate intermediate-metallicity stars. 
Main Sequence (MS) dwarfs were excluded both by a sloped cut at fainter magnitudes, parallel 
to the MS, which excludes all but the nearest (d<63 pc) dwarfs, and by co nservative kinemat ical 
criteria imposing a stellar velocity lower than the local escape velocity (see lGirard et al.ll2006l . for 
more details). 

We collected high-resolution Echelle spectra for two-thirds of the Girard et al.'s sample (824 
stars). All the stars in the brighter half of the sample were spectroscopically observed, while the 
completeness decreases with magnitude for fainter objects. We found a high contamination (~35%) 
by nearby dwarf in the faintest end of the distribution, hence stars with Kg > 9.5 were given lower 
priority after the first observations. The distribution of proper motions and colors of the observed 
sub-sample is shown in Figure [H The comparison with the complete sample reveals that no selection 
effect was introduced: the slight overabundance of red stars in the observed sample is due to the 
higher completeness at brighter magnitudes, where stars are on average redder. 



Table 1: Log of the spectroscopic observations 



Run 


Date 




N. nights 


Instrument 


N. stars 


Coraliel 


2005, 


September 


4 


Coralie@Euler 


39 


FEROSl 


2005, 


September 


6 


FEROS@2.2 


183 


duPontl 


2005, 


October 


6 


Echelle@duPont 


168 


Coralie2 


2005, 


October 


3 


Coralie@Euler 


25 


FER0S2 


2006, 


August 


6 


FEROS@2.2 


161 


duPont2 


2006, 


September 


6 


Echelle@duPont 


119 


MIKEl 


2006, 


Aug.-Nov. 


5x0.5 


MIKE@Magellan 


77 


MIKE2 


2007, 


October 


2 


MIKE@Magellan 


52 
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As discussed later, we will restrict our investigation to stars with Galactic height \Z\ >1.3 kpc, 
to avoid a strong thin disk contamination. This reduced the sample under study to 514 stars. We 
visually inspected all the spectra, identifying 46 probable dwarf stars. As expected, they show on 
average low radial velocities but high U and V components, and only one is found at <9. We 
also found 22 stars with [Fe/H]< —1.5, as deduced from comparisons with synthetic spectra. They 
were considered probable halo contaminants and were excluded from furth er analysis. We note that 
the m etal-poor tail of the thick disk extends to much lower metallicities (jBeers fc Sommer-Larsen 



1995). but below this limi t thick disk stars are outnumbered by halo members by a factor of nearly 



ten (jChiba &: Beera 120001 ). Finally, we also excluded from analysis 34 stars that, in the Galactic 
cylindrical coordinate system, had velocity components (defined in Section 14. 2p outside the range 
|U| <300 km s-\ -500 <V<300 km s-\ and |W| <150 km . The cut in W was aimed to 
remove the residual halo contamination, while the other components were restricted to exclude 
stars with wrong distance or bad proper motion. The cut in V was offset toward negative values to 
take into account the vertical shear (see Section [5.ip . Our final sample thus comprised 412 stars. 



2.1. Distances 



The interstellar reddening E(i? — V) was derived for each target from the lSchlegel et al.l (119981 ) 
m aps, and the JC s magn itude and {J — Kg) color were de-reddened by means of the transformations 



of ICardelli et al.l (|l989l ). The distance to each star was then estimated through a color-absolute 
magnitude relation calibrated on 47 Tucanae, a disk globular cluster (|Zinnl Il985ll whose stellar 
popu lation is very similar to the Galactic thick disk both in age and metallicity (jWvse &: Gilmore 



20051 ). The fit of the cluster red giant branch yields the relation (L. Miller 2008, priv. comm.): 



Kg = -7.886 • (J - Kg) + 16.302, 



(1) 



which were translated into absolute magnitude and de-reddened color assuming (m— M) v=13.50ifc0.08 
and E(i?— y)=0.024zb0.004 for the cluster distance modulus and reddening, respectively (iGratton et al, 
2OO3I ). 



Inspecting the 2MASS photometric data of 47Tuc used to derive Equation ([T]), we found that 
the uncertainty on Mk is of the order of ~0.2 magnitudes. This is only a marginal contribution to 
the final error in distance, because the main source of uncertainty is the relatively wide range of 
age and metallicity covered by thick disk stars. Indeed, the 2MASS photometric errors have only 
negligible impact, because they do not exceed 0.03 magnitudes for our fainter targets {Kg ~10.6). 
We estimated the effect of the age an d metallicity di stribution on the derived absolute magnitudes 



by means of Yale-Yonsei isochrones (|Yi et al.ll2003l ). W e assumed a scatter of 0.3 dex in metal 



scarc ely-populated low- and high- m etallicity tails (IBeers &: Sommer-Larsen 
2OO7I ). and a scatter of 2 Gyr in age (|Bensbv et ahlbood : iFeltzing et aPbood : 



1995 




Bensbv et al. 


Reddv et al. 


2006 


)■ 



They were considered uncorrelated, because the age-metallicity relation for the thick disk is very 
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weak (jBensbv et al.l 120071 ) . We finally estimated the error on distance to be 
summing all the relevant contributions. 



--20%, quadratically 



Thin disk stars do not follow the age and metallicity distribution assumed to derive Equa- 
tion ([1]), and their distances should be systematically biased. Indeed, younger, more metal-rich red 
giants are intrinsically fainter than our estimate. Comparing the absolute magnitudes calculated 
from Equation ([T|) with Yale- Yonsei isochrone s bet ween 2 and 8 Gyr , and metallicity following 
the age-met allicity relation of iHavwoodl (|200ll ) and iRocha-Pinto et al.l (|2000l ). we found that the 
distance of thin disk stars would be overestimated by 10-20%. This systematic error is small, but 
it has a relevant consequence on the contamination of the sample, which is artificially stretched to 
larger heights from the plane. 



2.2. Thin disk contamination 

Our sample is contaminated by the thin disk, because the targets were selected through a 
color cut that efficiently eliminates low-metallicity stars, but which excludes only a small fraction 



of metal-rich objects, as it can be deduced from Figure 1 of iGirard et al.l ()2006l ) . In Figure [2] we 



plot an estimate of the fraction of thin disk stars in the sample, assuming 300 and 90 pc as thin 



and thick disk scale height, respectively, and a thick disk local normalization of 12% (|Juric et al 



20081 ). We also calculated the contamination in the case of a distance overestimate of 15% for thin 
disk objects. The curves can be considered an upper limit, because the color cut, unaccounted for 
in the calculation, could have removed part of the contamination. In the same figure we indicate 
the results of a rough empirical estimate, obtained fitting a double Gaussian to the distribution of 
the vertical velocity component (W) in five bins of 50 stars. The dispersions of the two Gaussian 
components were kept fixed, and the only free parameter of the fit was the number of thin disk stars 
in the bin. The error bars show the results when varying aw in the range 14-16 and 35-40 km s^^ for 
the thin and the thick disk, respectively. This test was performed only for \Z\ ~l-2 kpc, where both 
populations contribute by more than 20% of the objects, i.e. 10 stars in each bin. The empirical 
errorbars agree well with the curve obtained when the distance bias is considered, although the 
observed contamination fall-off with \Z\ seems steeper than the modeled one. 

In conclusion, our sample is contaminated by thin disk stars, whose distances are overestimated. 
This affects even the estimate of their kinematics, which is also biased. Therefore, we will not 
analyze the results for the thin disk, and we will limit our study to \Z\ >1.3 kpc, to avoid the 
strong contamination (>50%) of the nearest subsample. 



2.3. Halo contamination 



Girard et al.l (|2006l ) estimated that the halo contamination in the sample should be about 



70, 



hence we would expect 41 halo members among the 514 targets with \Z\ >1.3 kpc. The cut in W 
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Fig. 1. — Left panel: color distribution of the spectroscopically observed stars (thick line) and of 

the complete sample (shaded histogram); Right panel: proper motion distribution of the spectro- 
scopically observed stars (black dots) and of the complete sample (grey dots) . 




Fig. 2. — Estimated fraction of thin disk stars in the sample. The curves indicate the expectation 
of the model described in the text, when a distance overestimate of 15% for thin disk is considered 
(dashed curve) or not (thick curve). The errorbars show the results of our empirical estimate. 
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removed 10 probable halo contaminants and, according to the statistics of IChiba &: Beerd (|200Cll ). 
18 of the excluded low-metalhcity targets are expected to be halo stars. The cuts in U and V could 
also have removed some halo objects, but the effect of this cut is harder to quantify, because an 
unknown fraction of the outliers are probably objects with problematic proper motion or wrong 
distance. The residual halo contamination in our sample of 412 targets should therefore be lower 
than ~3% (13 stars), and it can be safely neglected. 



3. OBSERVATIONS AND DATA REDUCTION 

The spectra were collected during six observing runs between 2005 and 2007, at La Silla and 
Las Campanas ob servatories. The details of the observations and data reduction were presented in 



Moni BidinI (|2009l ). and we will outline here only the most relevant points. Four different telescopes 
were used, because the stars span a wide range in magnitude {V=5-16), but the instruments had 
similar characteristics. The basic information of the runs is given in Table [TJ The spectra covered 
the whole visual range from the atmospheric cutoff on the blue side to ~9000 A, except for Coralie 
data which only reached 6800 A on the red edge. We did not reduce the spectra collected with 
the MIKE blue arm because of too low stellar flux, and MIKE data were thus limited to 4850 A 
on the blue end. The spectral resolution varied between 32 000 and 50 000, depending on the 
spectrograph. In particular, the 0'.'7 slit was always used for MIKE runs (R=32 000), while at 
duPont both the ^.'75 and 1" slits were used (R=40 000 and 30 000, respectively), depending on 
weather conditions. During each run we secured the spectra of up to seventeen bright stars with 
accurate parameters from the literature (radial and rotational velocities, temperature, gravity, 
metallicity) . Exposure times were chosen according to weather conditions, to reach S/N=70-100 
for all the targets. Comparison lamp arcs were acquired only during daytime calibrations for the 
fiber spectrographs Coralie and FEROS. The second fiber of these two instruments was allocated 
to the sky background, because the spectra were usually collected next to full moon. Lamp arcs 
were collected at intervals of about 2 hours during duPont and MIKE runs, and each spectrum 
was calibrated with the average of the two lamps temporally enclosing it, linearly weighted by the 
temporal distance between science and calibration frames. 

Spectra were reduced with standard IRAllfl tasks, and we kept the procedures as much as 
possible identical for all the data. We a posteriori verified that the reduced FEROS and Coralie 
spectra were of the same quality as the products of their dedicated pipelines. We analyzed the dark 
frames collected for all the runs, but we always found a negligible instrumental dark current and no 
clear 2D pattern, hence no dark correction was applied. The wavelength calibration lamp spectra 
were extracted in the same position on the CCD as science targets, to avoid systematics introduced 



■^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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by the curvature of the orders. The spectrum of a bright, well-exposed star was used to trace the 
orders on the CCD in each observing night, allowing only for rigid s hifts among the frames. Then, 
the spectra were extracted with an optimum extraction algorithm (|Hornelll986l ). normalized, and 
finally merged. 



MEASUREMENTS 



4.1. Radial velocities 



T he radial velocity (RV ) for all the targets was measured with a cross-correlation (CC) tech- 
nique (|Tonrv^^avi] 1 197d}_ as implemented in the IRAF fxcor task. The procedure was described 



in detail in iMoni BidinI (|2009l ). and we give only a brief summary here. The spectrum of each 
star was cross-correlated with three template stars observed in the same run, namely HD 180540, 
HD 223559, and HD 213893 (this last replaced by HD 224834 for Coralie data), encompassing the 
temperature range of the targ ets. Their RV was fixed by a CC with a synthetic spectrum generated 
with the SPECTRU1V0 code ( Grav Sz Corballv 1994), because we found poor agreement between 
the available literature sources. The analysis of the seventeen standard stars acquired during obser- 
vations, and of the solar spectra collected each night, revealed that the RV zero-point thus defined 
was biased by 0.3-0.7 km s~^, depending on the instrument, and this offset was removed. 

The three measurements were averaged, although they never differed by more than 0.2 km s~^. 
The CC was restricted to the range 4600-6800 A (5000-6800 A for MIKE data), to avoid the poor- 
quality blue section and the telluric bands at longer wavelengths. All RVs were corrected to 
heliocentric velocities, then the position of the telluric band at 6875 A was used to correct the 
RVs for instrumental effects, mainly caused by a displacement on CCD between the lamp and 
sci ence spectra, and an of f-center position of the star inside the slit (see, for example, the analysis 
of iMoni Bidin et al.ll2006l ). Corrections of up to 2 km s~^ were applied, but with little scatter 
(~0.5 km s~^) within each observing night. 

The final RV errors were estimated as the quadratic sum of the five relevant uncertainties 
introduced in the procedure: the CC and wavelength calibration error, the uncertainty on the zero- 
point and its offset, and the error on the estimate of the correction for instrumental effects. The 
resulting errors were typically in the range 0.4-0.7 km s"^ The final RVs of all the 824 stars wih 
be made available at the CDS websit^. 

We found 211 previous RV measurements for the stars in our sample. Our results agree excel- 
lently with literature sources: the mean difference (in the sense ours— literature) is 0.4ib2.7 km s~^, 
where the uncertainty is the rms of the differences. The comparison with the works that share 



^http: / /www.phys. appstate.edu/spectrum/spectrum. html 
^http://cdsweb. u-strasbg.fr/ 
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with us more than five stars in com mon is given in Table El The mean RV difference is always 
of the order of 1 km s~^, except for iMaiewski et alj ( 20041 ). whose RVs are higher than ours by 
a non-negligible amount (~5 km s~^). However, the number of stars in common is too small to 
conclude that this offset is significant. 



4.2. Galactic velocities 



The proper motion, radial velocity, and distance of each target were combined to derive its 
({/, V, W) velocity components in the Galactic cylindrical reference frame, where U is positive 
toward the Galactic center, V is directed in the sense of Galactic rotation, and W points toward 
the North Galactic Pole. The error on these velocities was derived propagating the uncertainty on 
the observed quantities. We assumed a proper motion error of 3 mas yr~^ for all the stars, as this 
value is more re alistic than the un certainties quoted in the SPM3 catalog (T. Girard, 2009, priv. 
comm. ; see also iGirard et al.l l2006l ) . The veloc ities were corrected subt ra cting the solar peculia r 



motion {Uq,Vq,Wq)={W.O, 5.1, 7.2) km s'^ (IPehnen fc Binnev 



1998 



Schonrich et all pOld ) 



recently proposed the revised values {Uq,Vq,Wq)={11.0, 12.2, 7.3) km s~"^, but we preferred to 
adopt the older ones, for sake of continuity with previous works. In any case, the definition of the 
solar motion does not affect the velocity dispersions nor the off-diagonal terms of the dispersion 
matrix (Equation [2]) , while the effects on the mean values are discussed in Section 15.11 

The sample was then divided into several bins, in which we calculated the mean velocities, the 
dispersions {a\j , ay , a^j^) , and the non-diagonal terms of the dispersion matrix 

1 



(iV-1) 



(2) 



where the sum is extended to all the stars in the bin, and a, (3 = U, V, W. The results are given 
in Table El The bins were defined by the requirement that their centers were equally spaced by 
0.1 kpc from |Z|=1.5 to 4.5 kpc, to uniformly sample the variation with \Z\ of the kinematical 



Table 2: Mean RV difference between this work and the quoted reference for the stars in common. 



Reference 


N. stars 


Arv 






km s~^ 


Flvnn & Freeman fl993) 


145 


0.8±2.5 


Zwitter et al. f 20081 


9 


-1.2±1.5 


Beers &: Sommer-Larsen f 19951 


9 


0.4±1.8 


Maiewski et al. (2004) 


8 


-4.7±2.1 


Duflot et al. f 19951 


6 


-1.1±3.8 


Jones (1972) 


6 


-0.7±0.9 


Evans ('19701 


6 


1.8±4.2 
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quantities, while the width was imposed by the fixed number of stars per bin (see below). This 
implied a partial overlap of the bin s, increasing with distari ce from the plane due to the decreasing 



density of observed stars. However, iMoni Bidin et al.l ([20101) have shown that the binning definition 



does not alter the results, and in fact our results are very similar to that work, despite the very 
different binning schemes. The bins with \Z\ >2.4 kpc comprised 50 stars each, while at lower 
heights, where the number of observed stars is larger, the bin size was increased to 80 targets for 
2.1< \Z\ <2.4 kpc, and 100 targets for \Z\ <2.1 kpc. We thus ensured a good statistic of thick disk 
stars in bins contaminated by the thin disk. 

In each bin, the mean value and dispersion of each velocity component were measured b y 



means of the analysis of the corresponding probability plot (|Hamakerlll978l : iLutz &: Hansonlll992l ). 
This powerful tool is highly insensitive to outliers, and it can be employed even in case of poorly 
populated bins. The data are first ordered with increasing value of the velocity. Each point i is 
then assigned a value cr, equal to the deviation from the mean expected for the i*'^ point of a normal 
distribution, in units of the standard deviation. When the underlying distribution is Gaussian, the 
data follow a linear relation in the a-velocity plane, whose intercept and slope are given by the 
mean value and the dispersion, respectively. The probability plots were fitted with a straight line, 
weighting each data point by its error, two examples are shown in Figure [3l The intrinsic velocity 
dispersions were then obtained quadratically subtracting the mean velocity error in the bin. The 
wings of the distribution showing any deviation from linearity were excluded from the fit, suspected 
to be affected by objects with wrong distance or problematic proper motion, or by residual halo 
members. 

The analysis of artificial data samples, generated through Monte-Carlo simulations and an- 
alyzed as the real data, indicated that the velocity dispersion is underestimated by 1-2 km s^^ 
when the thin disk contamination approaches 10%. The contamination was therefore neglected for 
\Z\ >2.5 kpc, where it is e xpect ed to be lower than 5% (i.e. two stars per bin, see Figure [2]). 



Following iBochanski et al.l (|2007l ). in the presence of a non- negligible thin disk contamination 
(\Z\ <2.5 kpc) we derived the mean velocity and the dispersion of thick disk stars fitting only 
the wings of the probability plot, and only the negative one for V . The cut excluded the ranges 
|U| <50 km s~^, |W| <30 km s~^, and V > —60 km s~^, approximatively 1.5 times the expected 
thin disk dispersion, thus removing about 90% of the contaminants. The two wings of the proba- 
bility plot should have the same slope but a different intercept, and the mean velocity was obtained 
averaging the two intercepts obtained from the fit. Although the mean values thus derived agree 
well with the trends observed in the contamination- free bins at higher |Z|, we regard them as less 
reliable. 

The formal errors of the least-square procedure used to obtain the intercept and the slope in 
the probability plots mainly indicate the deviations from Gaussianity of the observed distribution. 
Hence, they are not a good estimate of the real uncertainties. The errors were therefore derived by 
means of Monte-Carlo simulations. For each bin, we generated one thousand artificial samples of 50 
stars, changing the proper motion, distance, and RV of each star assuming a Gaussian random noise 
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Fig. 3. — Examples of probability plots used to derive the mean value and the dispersion of the 
velocity components. The velocity is plotted against the cr-value, assigned to each datum as de- 
scribed in the text. The mean value and the dispersion are given, respectively, by the intercept and 
the slope of the linear fit to the points in the plot. Upper panel: bin of 50 stars with \Z\ >2.5 kpc. 
The line shows the fit, and the empty dots indicate the data excluded from the fit. Lower panel: 
bin of 100 stars in the region contaminated by thin disk stars (\Z\ <2.5 kpc). Only the wings of 
the distribution were fitted, and the data not used in the procedure are shown as small dots. 
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with dispersion equal to the estimated errors (3 mas yr~^, 20%, and 0.4-0.7 km s~^, respectively). 
The measurements were repeated in each artificial bin as done for the real data, and the dispersion 
of these artificial measurements was assumed as an estimate of the observational errors. The final 
errors are given in Table [3l 

The non-diagonal terms of the dispersion matrix, cTuy, 0"^^, and cTy^, were calculated by 
means of Equation [2j The errors were estimated from Monte-Carlo simulations, as done for the 
velocity dispersions. For each bin, one thousand artificial data sets of 50 values were simulated, 
drawn from a distribution characterized by the same dispersion matrix as the real data. In each 
simulation, a Gaussian random noise with a dispersion equal to the observational errors was added, 
and the non-diagonal terms were estimated by means of Equation [2l The errors were defined by 
the scatter of these simulated measurements. With the same simulations we also evaluated the 



systematic error introduced by the observational uncertainties (e.g. ISiebert et al.ll2008l ). comparing 
the results when the random noise was considered or not. We found that the expected systematic 
thus introduced is negligible, being about one order of magnitude lower than the uncertainties on the 
measured values. Therefore, we did not correct the observed non-diagonal terms for observational 
errors. 



RESULTS 



5.1. Mean velocities 



The trend of the mean velocities with distance from the plane is shown in Figure [H W 
is consistent with no significant departure from zero in the whole range. Some deviations are 
observed for \Z\ <2.5 kpc, but in this range the results are affected by large errors and, as already 
commented, not very trustworthy. This result agrees w ith previous studies, that found no net 
vertical motion up to variou s kpc from the Galactic plane (jBochanski et al.ll2007l : iBond et al.ll20ld : 



Casetti-Dinescu et al 



2011 



). On the contrary, U{Z) has a more complex behavior: between |Z|=1.5 
and 3 kpc we detect a non-null mean velocity directed toward the Galactic anticenter, with an 
average value C/=— 19±3 km s^^. Beyond 3 kpc from the plane U abruptly increases and changes 
sign, and the net velocity between 3.5 and 4.5 kpc is 15ib2 km s^^ toward the Galactic center. 
Another possible interpretation of the results is that U linearly increases with \Z\, and a S-shaped 
feature is superimposed to this trend between 2.5 and 3.5 kpc, as observed for a\j and other 
kinematical quantities (see Section [5. 5p . After the exclusion o f this range, the f it retu rns a vertical 
increase of 15=b2 km s^^ kpc^^. Adopting the new values of ISchonrich et al.l ( 2010l ) for the solar 
peculiar motion causes not a relevant change on the results, as U would be higher by about 1 km s^^ 
only. 

Previous studies have shown evidence that stars more distant than 1 kpc from th e plane show 
a non-null net radial motion of the order of ~9 km toward the Galactic anticenter ( Smith et al. 

Our results between |Z|=1.5 and 3 kpc 



2009 



Rvbka Sz Yatsenko 



201C 



Casetti-Dinescu et al. 



2011 
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Table 3: Measured kinematical quantities and associated uncertainties in each bin. 



\z\ 


U 


V 


W 




cry 


0"W 


2 


9 


2 


kpc 


km s 


km s 


km s 


km s 


km s 


km s 


km^ s 


1 9 9 

km^ s 


km^ s 


1.5 


-31±4 


-46±4 


17±1 


81±5 


57±5 


38±1 


570±310 


400±470 


10±210 


1.6 


-20±4 


-55±4 


5±1 


77±5 


59±5 


37±1 


840±290 


-180±460 


-20±220 


1.7 


-28±4 


-61±5 


3±1 


79±5 


64±5 


38±1 


930±310 


-130±510 


-40±250 


1.8 


-24±4 


-69±5 


8±1 


79±5 


61±6 


40±1 


810±330 


-760±500 


-80±240 


1.9 


-16±4 


-77±5 


3±1 


78±5 


60±5 


40±1 


850±330 


-1060±500 


-120±240 


2.0 


-10±5 


-85±4 


4±1 


83±6 


55±6 


39±1 


940±340 


30±460 


-90±210 


2.1 


-22d=5 


-78±5 


3±1 


80±6 


58±6 


38±1 


1060±320 


470±470 


-130±220 


2.2 


-16±5 


-76±7 


-13±1 


77±6 


63±7 


42±1 


690±330 


0±550 


230±260 


2.3 


-1±6 


-81±5 


-9±1 


81±7 


58±6 


40±1 


750±340 


-500±470 


10±230 


2.4 


-12±5 


-85±6 


-12±1 


80±6 


59±6 


40±1 


410±330 


-710±480 


180±240 


2.5 


-18±6 


-85±6 


-2±1 


78±7 


63±7 


42±1 


630±330 


-1870±600 


-40±270 


2.6 


-25±6 


-90±6 


1±1 


71±7 


66±7 


42±1 


530±310 


-2390±660 


-460±290 


2.7 


-23±6 


-98±6 


7±1 


72±7 


62±7 


39±1 


870±290 


-1150±520 


-1360±300 


2.8 


-29±6 


-95±7 


-2±1 


76±7 


62±7 


41±1 


880±320 


-1740±580 


-1310±300 


2.9 


-17±6 


-115±7 


5±1 


83±7 


68±7 


40±1 


840±340 


-820±590 


-1230±300 


3.0 


-7±6 


-126±7 


8±1 


90±7 


67±8 


42±1 


1440±410 


-50±610 


-460±280 


3.1 


-6±7 


-129±8 


1±1 


101±8 


67±8 


43±1 


1610±450 


-1150±690 


-650±290 


3.2 


2±6 


-131±8 


-1±1 


99±8 


63±8 


42±1 


1470±440 


-1820±650 


-440±280 


3.3 


8±7 


-140±8 


Oil 


101±8 


66±8 


44±2 


2180±490 


-2610±740 


-310±300 


3.4 


12±7 


-140±8 


1±1 


98±8 


63±9 


43±2 


2100±480 


-1880±660 


260±280 


3.5 


18±7 


-135±9 


5±1 


95±9 


63±9 


44±2 


1950±460 


-3390±740 


260±280 


3.6 


29±7 


-137±8 


9±1 


101±9 


64±9 


44±2 


2290±510 


-2170±700 


360±290 


3.7 


7±7 


-132±9 


0±2 


91±9 


61±10 


44±2 


2180±460 


-2150±610 


610±280 


3.8 


15±7 


-134±10 


1±2 


92±9 


68±10 


43±2 


1990±450 


-2120±700 


-30±300 


3.9 


5±8 


-139±9 


-6±2 


94±10 


66±10 


43±2 


1540±430 


-520±620 


210±280 


4.0 


7±8 


-143±10 


-1±2 


93±9 


66±11 


41±2 


1910±420 


-100±610 


10±280 


4.1 


17±8 


-145±10 


1±2 


92±10 


68±10 


44±2 


1570±440 


720±640 


-30±300 


4.2 


22±8 


-142±11 


2±2 


95±10 


72±11 


45±2 


1480±460 


880±710 


10±320 


4.3 


12±9 


-148±10 


-2±2 


94±10 


69±11 


48±2 


2320±510 


1170±680 


180±330 


4.4 


17±9 


-152±11 


0±2 


96±10 


72±11 


46±2 


2330±500 


1370±710 


-50±330 


4.5 


19±10 


-158±11 


-2±2 


93±11 


76±11 


46±2 


2180±480 


3960±930 


240±360 
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Fig. 4. — Trend of mean velocity components (from top to bottom: radial, rotational, and vertical 
component) as a function of distance from the Galactic plane. Empty dots are used for the bins 
contaminated by the thin disk, where the me a surein ents ar e less reliable. Result s of p revious 
investigations are als o indi cated: 



Chiba fc Beer j pnod . COO), ISpagna et al.l poiol . SIO) 



Girard et al 



(|2006l . G06), ICasetti-Dinescu et all (|201ll . Dll 
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agree with their conclusion, although our mean value is larger by about a factor of two. The sign 
flip observed at lar ger \Z\, on the cont rary, has never been repor ted in the literature. Nevertheless, 



the plane, a nd 



Smith et al.l (120091) analyze halo stars . 



Casetti-Dinescu et al. ( 



2011 



Rvbka &: Yatsenkd (j2010l ) do not reach these heights above 



have too few stars in this range. Very interestingly, 
Bond et al.l (|2010l ) detect a small positive mean radial component out to 1^1 3 kpc, and a negative 
value in their more distant bins, alt hough the authors coii c lude that these results are consistent, 
within errors, with a null net motion. ICasetti-Dinescu et al.l (|201ll ) proposed that an inward motion 
of the Local Standard of Rest (LSR) of the order of 10 km s~^ is responsible for the observed non- 
null mean value of U, as suggested by the fact that this is observed among both disk and halo stars 
more distant than 1 kpc from the sun. This hypothesis cannot account for any vertical trend of U 
other than a constant non-zero value at any \Z\. Thus, while it is not contradicted by our results, 
at least another effect must be invoked. For example, an outward motion of ~25 km s~^ of the 
stars between 3.5 and 4.5 kpc, coupled with the mentioned LSR motion, could explain our results. 
This clumpy kinematical distribution would not be surprising, because perturbations produced by 
the bar an d the spiral arms a re expected to form groups of stars with velocity offset as large as 
50 km s"^ (jOuillen et al.lboiol ). 



Girard et al.l (|2006l ) studied the proper motions of our same stellar sample, and their results for 
U (Z) are overplotted to ours in the upper panel of Figure HI after correcting for the solar peculiar 
motion and changing the sign of U to account for the different direction of the reference axis. 
The agreement with our results is poor: while their nearly flat profile roughly coincides with our 
mean value (—7.7 km s~^), they did not detect any change of sign, nor a steep positive gradient. 
The different approach to the same data must have ca used the different r esults. For example 
the features observed by us could have been masked in I Girard et al.l (120061 ) by their smoothed, 
statistical distance determination, or by dwarf stars and halo contaminants, removed in our work. 

In the middle panel of Figure H] we show the measured trend of V{Z), compared to other 
results from the literature. The mean rotational velocity of thick disk st ars decreas e s wit h distance 
from the plane, due to the well-known vertical shear first detected by iMajewskil (119921 ). The fit 
of the data, after the exclusion of the less reliable bins at \Z\ <2.5 kpc, yields V{Z) = (—25 it 
12) - (30 ifc 3) ■ \Z\ (km s~^) . As can be seen in Figure HJ this solution e x cellen tly matches the 



results of iGirard et al.l (j2006l ). and even the data points of IChiba &: BeersI (j2000l ). at z < 2 kpc, 
are well aligned with the derived linear relation. Both these in v estiga tions measure a vertica l 
shear of —30 km kpc~^, as also recently found bv lBrown et al.l (j2008l ). and lBond et al.l ( 2010l ). 
Had we fitted all the data points down to |Z|=1.5 kpc, we would have found a steeper slope 
(-35.1 ± 1.8 km s "^ kpc~M, at the up per li mit of the range sp anned by literature values, which 



vary from — 16ib4 (jAllende Prieto et al 



20061 ) to -36zfcl km s~^ (ICarollo et al.ll20ld ). The revised 



values for the solar peculiar motion proposed by ISchonrich et al.l (|2010l ) cause an upward revision 
of the results by 7.1 km s~^. 



Casetti-Dinescu et al.l (|201ll ) and lSpagna et al.l (j20ld ) measured a shallower slope (— 25.2ib2.1 



and — 19ib2 km s ^ kpc ^, respectively) between 0.7 and 2.8 kpc. While our results are compat- 



- 16 - 







-50 



A 
> 
V 



- 100 



150 







^ coo 
■ sio 

^ DU 

^ this work 




J \ I \ I L 



5 



Fig. 5. — Rotational velocity of thick disk s tars as a function of distance from the plane. Full dots: 
our work; asterisks: IChiba. &: Beerd (j200d . COO); squares: ISpagna et al.l (j2010l . SIO); triangles: 



Casetti-Dinescu et al 



(2011 



, Dll). The dotted curve indicates the power-law solution proposed by 
Bond et al.l (j2010l ). the thick curve is the analogous relation obtained from the fit of the plotted 
data points. 
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ible with theirs in our nearest bins, we find a lower mean rotational velocity beyond \Z\=1 kpc. 
Nevertheless, the different vertical rotational gradient found by these studies is not necessarily 
a disagreement, because they sample a different jZl-raii ge, and the underlying shear is not re- 
quired to be strictly linear. In fact, llvezic et alj (|2008l ) have proposed the non-linear relation 
VfZ) = —20.1 - 19.2 • 1-25 g-i from t he analysis of SPSS data. Th e combined data points of 
Chiba Beera (j200Cl ). ISpagna et al. ICasetti-Dinescu et al.l (120111 ) , and of the present work, 



closely follow this equation (see Figure [5]), and a fit in the form V i^Z^ = a ^ ■\Z\^ km s""*^ returns 
a very similar solution, with a = —22.5 km s~^, (3 = —22.2 km s^-*^ kpc^-*^, and 7 =1.23. Very 
noticeably, the results of four surveys finding a different linear relation are all well reproduced by 
a single non-linear curve proposed by an independent work. The underlying vertical trend of the 
thick disk rotational velocity is therefore most probably close to but not exactly linear, and better 
represented by a low-order power law. 



5.2. Velocity dispersions 



The vertical profile of the velocity dispersions is shown in F igure [6l coinpared to other works 
from the literature. We excluded from comparison the results of I Bond et al.l (|201Cl ). because they 
do not distinguish between different disk components, thus finding steeper gradients as a result 
of the mix of thin and thick disk stars. Caution must also be taken when comparing the results 
obtained in different ranges of |Z|, because the underlying trend of the dispersi ons is not necessarily 
linear, and the gradient can assume different values. For example, the models of lGirard et al.l (|2006l ) 
indicate that the vertical profile should be progressively flatter at increasing \Z\. 



The linear fit of the trends with \Z\ yields: 



av{Z) = (82.9 ± 3.2) + (6.3 ± 1.1) • {\Z\ - 2.5)km s" 



(3) 



(jv(Z) = (62.2 ± 3.1) + (4.1 ± 1.0) • (i^l - 2.5)km s"\ (4) 

aw{Z) = (40.6 ± 0.8) + (2.7 ± 0.3) • {\Z\ - 2.5)km s"\ (5) 

where \Z\ is in kpc. The quoted errors were obtained restricting the fit to a subset of eight non- 
overlapping bins, to avoid the unde restimate induced by th e correlation between the data points. 
The results are identical to those of iMoni Bidin et al.l (|2010l ). despite the different binning scheme, 
except for cry for which we derive a slightly smaller gradient, but the difference is not significant 
(0.4 km s~^ kpc~^). The gradients are small, and the three dispersions increase by only ~7% 
between 2.5 and 3.5 kpc. This explains why the change of the thick disk kinematics with distance 
from the plane has not been clearly detected for about two decades after its discovery. 



The results for ctw are the most precise and reliable, because ~90% of the information on W 
comes from RVs, whose errors are an order of magnitude smaller than those in proper motions. 
a\\i{Z) shows a small but clear and constant increase, with small scatter around the der i ved lin ear 
relation. The vertical gradient is smaller than the one found by ICasetti-Dinescu et al.l (|201ll ) by 
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Fig. 6. — Velocity dispersion as a function of distance from the plane (from top to bottom: radial, 
rotational, and vertical velocity dispersion). The thick curve indicates the linear fit of the data 
given in Equations ([3]) to ([5]). The empty dots, dashed and dotted lines are as Figure HI 
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about a factor of two, but the results are consistent at the la level. lYoss et al.l (|l987l ) propose an 
even steeper profile (~10 km s~^ kpc~^) in the range \Z\=0-2 kpc, but the same authors suspect 
that this could be due to contamination by halo stars increasing with Galactic height. 

The dispersion of the rotational velocity component, ay, also shows a clear increase with \Z\, 
although the data points are affec ted by larger errors and are m ore scattered. The vert ical gradient 
is smaller than the results of both lCasetti-Dinescu et al.l (|201ll ) and lGirard et al.l (|2006l ). by a factor 
of three and two, respectively, and the difference is at the 2a level in both cases. 

The vertical profile of a\j shows a very peculiar behavior, with large deviations from linearity 
between 2.5 and 3.5 kpc. This feature will be discussed in Section [5.51 Outside this range, the data 
points closely follow a linear relation, and the solution given in Equation [3] was obtained excluding 
this interval from the fit. The derived vertical gradient of au is very similar to that found by 
Girard et al.l (j2006l ). but their solution is offset tow ard higher values by about 10 km s~^, while the 
gradient measured bv lCasetti-Dinescu et al.l (|201ll ) is 4o" times higher than ours. 



In conclusion, we confirm that the velocity dispersions of the Galactic thick disk increase with 



distance from the plane, as suggested by previous investigations (lYoss et al 



19971 : ISoubiran et al.l 120031 : iGirard et al.l 120061 : llvezic et al.l 120081 : ICasetti-Dinescu et al 



( 


Yoss et al. 


1987: 


Yoss & Griffin 


Casetti-Dinescu et al. 


2011 


), but 



we derive vertical gradients that are in general smaller than other studies. The differences with 
Casetti-Dinescu et al.l (120111 ) can be at least in part due to the aforementioned expected decrease 
of the gradient with distance from the plane, because the results for aw are consistent, and the 
solutions propos ed for the other two components overlap in the Z-range in common (1.5-2.2 kpc). 
On the contrary, iGirard et al.l (|2006l ) studied our same sample, and the different data analysis must 
be the cause of the disagreement. Since they measured both higher dispersions and steeper vertical 
gradients than us, their results could have been affected by the dwarf stars and halo contaminators 
removed by us, and/or the thin disk contamination, that they did not take into account. 

The three dispersions increase with distance from the plane proportionally, and the anisotropy 
is approximatively constant at all \Z\. In Figure [7| the vertical profile of the ratios a\]/a^ and 
a\]/a\i is shown, where the errors come from propagation of the uncertainties on the dispersions. 

The mean value of uu/cw is 2.08ib0.06, where the error was estimated from the statistical 
error on the mean for eight uncorrelated bins, as before. The linear fit of the data points returns 
a negligible gradient (0.09±0.07 kpc^^), and t he data are consi stent with a fiat profile. Literatur e 
values for au/cw span a wide range from 1.1 (IVeltz et al.ll2008l ) to 1.9 (e.g.. IVallenari et al.ll2006l ). 
and our resu lts are at the upper end of th is distribution. As shown in the upper panel of Figure [3 
the results of lCasetti-Dinescu et al.l (|201ll ) agree with our measurements in the Z-range in common, 
but they deviate from our measurements if extrapolated to higher distance from the plane. 

The results for a\]/a\ also show no significant gradient (0.06ib0.05 kpc~^) and a mean value 
of 1.36ib0.05. T his is lower than the valu e predicted by the epicyclic approximation (1.49, cfr. 
Equation 3-76 of lBinney &: Merrifieldlll998l ). indicating that the population under analysis cannot 
be assumed in nearly circular orbits. Our results are well within the range spanned by literature 
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7. — Ratio of dispersions as a function of distance from the plane. Upper panel: au/ay, lower 
1: (Tu/cw- The empty dots, dashed and dotted curves are as in Figure HI 



- 21 - 



value s, which vary from f«l (e.g.. lChiba &: Beersll200d:ICarono et al.ll20ld). t o «1.6 (jSoubiran et al 



20031 ). As in the case of au/a-w, the results of ICasetti-Dinescu et al.l (j201ll ) overlap with our data 
points in the range in common, but th eir extrapolation to higher \Z\ do not. The analytical 
expressions given bv iGirard et al.l (j2006l ) for a\j and ay, on the contrary, return a value of their 
ratio much lower than our data points at any Galactic height. This is probably due to their much 
higher ay , as shown in the middle panel of Figure El 



5.3. Orientation of the dispersion ellipsoid 



The non-diagonal term of the dispersion matrix gives a measurement of the correlation 
between the i-th and j-th velocity components. In fact, if two principal axis of the dispersion 
ellipsoid are aligned with the i-th and j-th axis, we have crf~0, and the two velocity components are 
uncorrelated. The orientation of the dispersio n ellipsoid of an old, dynamically re l axed population is 



related to the shape of the Galactic potential (jLvnden-Bell 



1962 



1991 


). Moreover, 


Bienavme ( 


2009) and 


Siebert et al. 



Ollongren 


1962 


Hori & Liu 


1963 




in the disk ( 


Kuiiken &: Tremaine 



ellipsoid allows to estimate the flattening of the dark halo. However, calculating the expectation 
of a specific halo model through integration of orbits is beyond the scope of the paper, and here 
we will only present the observational results, for use of future works. 

Our results for the three cross-terms are shown in Figure [8l The profile of a'^y closely follows a 
decreasing linear relation up to \Z\=A kpc, but abruptly jumps to positive values at larger heights. 
This is probably due to the sign fiip of U analyzed in Section [5.51 In any case, the results indicate 
that the U and V velocity components are correlated, and the velocity ellipsoid is titled in the 
radial-longitudinal plane. The rotation angle, i.e. the vertex deviation ip, was calculated in each 
bin by means of the relation: 



1 



—arctgl 



2a 



uv 



(6) 



(jAmendt &: Cuddefordlll99lh . The results are shown in the upper panel of Figure [9j We measure 
a non-negligible vertex deviation, increasing from nearly zero at |Z|=1.5 kpc to ~ 20° at 3.5 kpc. 
The linear fit in the range 1.5-4 kpc yields the relation ip = —1.0 + 5.4 • \Z\. Our results agree 
with ICasetti-Dinescu et al.l (j201ll ). who measured = 8.2 ± 3.2° at Z=l.l kpc. Previous investi- 
gations evidenced that the vertex d eviation decreases from 20° for young: populations to nearly 



zero for old , meta l-poor disk stars (iBienaymg Il999l : iDehnen &: Binney 



1998 



Soubiran et al 



2003 



Fuchs et al.l |2009|) . Our results are not at variance with this conclusion if the vertex deviation 
increases with \Z\, as suggested by our observations, because previous studies were limited to small 
Galactic heights. 

o^y^J^ shows no significant deviation from zero in the whole range of \Z\, except between 2.5 
and 3.5 kpc, as discussed in Section 15. 5i This indicates that the V and W velocity components 
are not correlated. On the contrary, ct'^t^ is significantly different from zero, and steadily increases 
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Fig. 8. — Vertical trend of the non-diagonal terms of the dispersion matrix, crfj^^, cr^y, and CTy-^^ 
(from top to bottom). The empty dots are used for bins contaminated by the thin disk. In 
the upper panel, the das hed line indicates the trend of the analytical expression proposed by 
Kuiiken k Gilmorel ^Wsd ). 
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with \Z\. Proposing a reliable linear express ion for (j^-^y^j{Z^ is not s traigh tforward, because it 
shows irregular fluctuations. As expected by iBinnev &: Spergell (|l983l ) and iBinnev &: Merrifield 



(|l998l ). cTuw is always bound between zero and (y^y^ ^^y^ = {U^ — W^) ■ (Z/R), the value assumed 



when the velocity ellipsoid is aligned with the spherical coordinate system, because we find that 
'^uw.max 4500 km^ s~^ at any Galactic height. We also find that the expression proposed by 



Kuijken &: Gilmord (|l989l ). obtained under the assumption that the ellipsoid points toward the 
Galactic center and its axis ratio is constant in the spherical coordinate system, is a relatively good 
approximation of the measured value up to about 3 kpc from the plane, but the predicted slope is 
too shallow and the agreement with observations degrades with \Z\ (see Figure [8]). 

The tilt angle in the U-W plane can be calculated, analogously to the vertex deviation, from 
the Equation: 



a 



1 

-—arctg 



2a 



uw 



u '^w 

and the results are shown in Figure [9l We measure a slight increase of a with \Z\, and the lin- 
ear fit of the data, after the exclusion of two deviating points, yields the relat ion a(|Z|) = 9?6 + 
2?4 • \(\ Z\/kpc) — 2) 1 . The results agree with pre vious works at [Z[=l-1.5 kpc (iSiebert et al.ll2008l : 



201C 



Casetti-Dinescu et al.l l201ll ) , but not with the measurements of 



Fuchs et al. 



Car olio et al 

(I2OO9I ) at \Z\ <1 kpc, whose sample is most probably dominated by thin disk stars. The dashed 



curve in Figure [9] indicates the value of a when the dispersion ellipsoid points toward the Galactic 
center, assuming R0=8 kpc. The tilt angle is constantly lower than this in the range 1.5-4.5 kpc, 
and the ellipsoid thus is directed toward a point located behind the Galactic center, at a Galacto- 
centric dist ance increasing with \Z\ from Ro=2.3 to 9 kpc. O r bit integration studies indeed predict 



this r esult (jKuiiken fc Gilmord 1 19891 : IBinnev Sz Spergel 



1983 



Kent Sz de Zeew 



1991 



Shapiro et al 



2OO3I ). but the inclination is noticeably higher than the expectations, because theoretical calcula- 
tions return Ro=5-10 kpc at \Z\=1.1 kpc. This could indicate that the Galactic potential used in 
these studie s needs to be refine d to match the observations. On the contrary, the increase of a{\Z\) 
modeled by iBond et al.l (I2OIOI ) is too steep, and their expectation a(3.5 kpc)=26° is at variance 
with our results. 



5.4. Modified Newtonian Dynamics 



Bienayme et al.l (120091 ) showed that the vertical trend of the tilt angle is an excellent ob- 



servational signature of the underlying gr avity law, and it can be used as a test for the Modi- 
fied Newtonian Dynamics theory (MOND. lMilgromlll983l ). because its expectation diverges from 
that of the New t onian dynamics with distance from the plane. Unfortunately, the calculations of 
Bienavme et al.l (|2009l ) are limited to \Z\ <2 kpc, and the ra nge of overlap with our data is very 
narrow. The value of a expected bv iBienavme et al.l (|2009l ) at \Z\=2 kpc is 12° for Newtonian 
dynamics and 10° for MOND. We find a{2 kpc)=9?5 it 3?6 and the mean of the five measurements 
in the range 1.8-2.2 kpc is 10?0 ± 0.5°, in excellent agreement with MOND expectations. The 
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Fig. 9. — Measured vertex deviation (upper panel) and tilt angle (lower panel) as a function of 
distance from the plane. The empty dots are used for the bins contaminated by thin disk. In 
the lower panel the dotted line shows the tilt angle of a dispersion ellipsoid aligned with the 
spherical coordinate system , at an y \Z\, and previous literature measurement s are also shown: 



Casetti-Dinescu et al.l (|201ll . Dlll. lCarollo et al.1 (|20ld . CIO), and lSiebert et al.1 (|2008l . SOS). 
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linear fit presented in Section [5.31 also implies a(2 kpc) =9?6. However, •we der ive a vertical gra- 
dient (2?4 kpc~^) much shallower than the predictions of lBienavme et alj ( 2003 ) for both MOND 
(5° kpc~^) and Newtonian dynamics (6 ° kpc~^). In conclusion, our observations agree better with 
the models that iBienavme et alJ (120091) der i ved fo r MOND than with those from Newtonian dy- 
namics, but the extension of lBienavme et al.1 (|2009l ) calculations to higher \Z\, and to specific initial 
conditions for the thick disk, are required to perform a reliable test of the gravitational law. 



5.5. Subtructures in the thick disk 

The vertical profile of cru{Z) shows a puzzling behavior between 2.5 and 3.5 kpc from the 
Galactic plane (top panel. Figure [6]), where the dispersion deviates from the linear trend. A 
similar feature is clearly observable in the trend of other quantities related to the radial velocity 



component U, like U and cr?f^ (as already noted bv lMoni Bidin et al.ll2010l ). but deviations from 



the linear trend in the range |Z|=2.5-3 kpc could be present even in other profiles, such as V{Z), 
W{Z), ay^{Z), and a'^^^{Z). The origin of this behavior is unclear. It is possible that a group 
of comoving stars, forming a sub-structure in the Galactic thick disk, is affecting the measured 
kinematics between 2.5 and 3.5 kpc from the plane. In particular, the similarity of the profile of 
two totally unrelated quantity such as Oyy^{Z) and (T\]{Z) is instructive, because it excludes the 
possibility that this behavior is only due to some bad measurements. In any case, the existence of 
a kinematical substructure among our stars cannot be claimed on the basis of these results only. 



5.6. Radial behavior of cju 

Our targets are distributed in a very narrow range of R, and the data do not provide any direct 
information about the radial behavior of the kinematical quantities. Nevertheless, some indirect 
indication can be derived. Manipulating the radial Jeans equation in cylindrical coordinates in 
steady state, with the radial component of the force expressed as Fr = — f^/R, we obtain: 



dal UW duw vl , 4 + y 2^1 1 . 



dR hz,p dZ R R ^Vi? Hr, 

2 

where we also assumed that U <C 0"^, d\J /dR=0, and that the density decays exponentially with 
both R and Z, with scale length hji^p and scale height Inserting the observed quantities 

in the right hand side of Equation [5] we can thus estimate, as an exercise, the radial behavior 
of fJu. The large uncertainties involved prevent a precise measurement, but the comparison with 
the expectations provides a consistency check for the kinematical data here presented, because an 
unphysical result could indicate a problem wi th them. The res ults are shown in Figure [TOl where 



we assumed hR^p=3.6 kpc and hz,p=0-9 kpc (jjuric et al.ll2008l ). R0=8 kpc, and Vc=220 km s ^ 



The data points are scattered around the mean value —10 it 1 km s ^ kpc ^, following the behavior 
of (Tu(Z) discussed in Section 15.51 with no clear vertical trend. In the nearest bins the results 
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Fig. 10. — Estimated radial derivative of a\j as a function of distance from the plane. The empty 
dots are used for the bins c ontaminated by thin disk. The squares correspond to the results of 



Casetti-Dinescu et al.l (|201ll ). The dotted line indicates the theoretical expectation assuming a 



radially constant anisotropy. 
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are partially consist ent, but more negative, than what found by ICasetti-Dinescu et alj (|201ll ) and 
Neese Yoss ( IQsd. —3.8 ± 0.6 km s~^ kpc~^), althou gh Neese & Yoss refer to the thin disk. The 



observations show (jvan der Kruit &: Searldll98ll . Il982l ') that aw should exponentially decay in the 
radial direction, with the same scale length of the mass de nsity. This is often con s idere d valid even 
for a"u, under the assumption of a constant anisotropy, and lCuddeford &: Amendtl (jl992l ) found that 
this should be the best approximation at the solar position. The results also roughly agree with the 
expectation of this model, as shown in Figure [101 In conclusion, the presented kinematics implies 
a radial behavior of a\j consistent with both previous observations and with the expectations of 
theoretical predictions. 



5.7. Comparison with models of thick disk formation 



The presence of kinematical gradients in the Galacti c thick disk is a powerful diagnostic to 
discriminate between the various models of its formation ([Majewskil Il993l ) . but unfortunately the 
model expectations concerning the vertical gradient have not been investigated so far. It would 
be interesting to compare, in the near future, the predictions of the various models with the 
observed trend with Z of the velocity dispersions. On the contrary, the radial gradient of the 
dispersions has been modele d in the context of thick disk formation through dis k heating by the 



merging of minor satellites (IVillalobos &: Helmil 120081 : iBekki &: Tsujimotd l201ll ) . The observa- 



tions of ICasetti-Dinescu et al.l (|201ll ) agree with the results of these simulations, provided that a 
low inclination orbit of the merging satellite is assumed. A similar conclusion is drawn compar- 
ing our measurements of cju/cw with the theoretical expectations of the merging scenario. In 
fact, the simulations showed that this ratio i s strongly linked to the inclination angle of the orbit 
(|Villalobos &: Helmill2008l : lPurcell et al.ll2009l). and our resu lt ((Tu/<7w=2.08) favors a small inclina- 
tion angle, i ~0-30° (see Figure 15 of lVillalobos et al.ll2010l ). A low inclination orbit of the infalling 
satellite is also required to reproduce the rat her large vertical gradient of the rotational velocity . 



found in this work as in previous investigations (IVillalobos &: Helmill2008l : lBekki &: Tsuiimotol[201ll 



On the contrary, the radial migration model predicts a much shallower gradient (—17 km s ^ kpc 



Loebman et al 



2011 



), incompatible with the observations. 



In conclusion, the models of thick disk formation have so far provided only fragmentary predic- 
tions about its kinematical properties, and our observations cannot still be used to fully discriminate 
between them. However, we find that our results are consistent with the scenario where the thick 
disk formed through dynamical heating of a pre-existing Galactic disk, induced by the merging of 
a minor satellite. Moreover, all the kinematical evidence shows that, if this is the correct model, a 
low-latitude (< 30°) merging event is strongly preferred. 
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6. CONCLUSIONS 

We have analyzed a sample of ~400 thick disk stars, measuring the variation of their kine- 
matical properties as a function of distance from the Galactic plane, from |Z|=1.5 to 4.5 kpc. Our 
results can be summarized as follows: 

• While the mean vertical velocity component W shows no significant deviation from zero in 
the whole range, between 1.5 and 3 kpc, we find a net radial motion of about 20 km s~^ 
directed toward the Galactic anticenter. Other authors have recently found evidence for a 
similar behavior, proposing an inward motion of the LSR, although our mean velocity is 
larger than their proposed value by a factor of two. However, we find that U changes sign 
for \Z\ >3 kpc, and that a radial motion of the LSR, although not excluded, cannot alone 
explain this behavior. 

• The mean rotational velocity of the thick disk decreases with distance from the Galactic plane, 
as found by many previous investigations. The linear fit of our data returns a gradient of 
—30 km s~^ kpc~^, our data points and the results of three other previous works are globally 
better repre sented by a power -law with index 1.25, very similar to what has recently been 



proposed by iBond et al.l (|20ld ) 



• All the velocity dispersions steadily increase with distance from the Galactic plane, closely 
following a linear relation. The gradients we found are, however, smaller than those proposed 
by previous works. 

• While the velocity dispersions increase with \Z\, the ratios ctu/<7w and au/ay show no sig- 
nificant vertical trend. The observations thus indicate a substantial constancy with |^| of the 
anisotropy. 

• We find a non-negligible vertex deviation, increasing with \Z\ from values close to zero to 
~ 20° at |Z|=3.5 kpc. This is consistent with previous investigations, which found a very 
small vertex deviation of old stellar population close to the Galactic plane. 

• The tilt angle steadily increases with distance from the Galactic plane. As expected, the 
orientation of the velocity ellipsoid in the U-W plane results, at any \Z\, intermediate between 
alignment with the cyl indrical and spherical coordinate systems. According to calculations by 



Bienavme et al.l (|2009l ). the tilt angle at \Z\=2 kpc coincides with the expectation of MOND, 
although the extension of their models to higher Galactic heights is required to perform a 
conclusive test of the underlying gravitational law. 

The vertical trend of many kinematical quantities show deviations from linearity between 2.5 
and 3.5 kpc. The origin of these features is unknown, but it could indicate the presence of a 
sub-structure at this Galactic height, such as a comoving group of stars. 
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• The results are fully consistent with the model of thick disk formation through dynamical 

heating of a pre-existing Galactic disk. If this is the correct scenario, a low inclination angle of 
the merging event is strongly preferred. However, not all the models proposed so far could be 
tested by our observations, and more simulations are required to obtain a detailed comparison 
able to discriminate between them. 
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